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The lower-critical dimension for the existence of the Ising spin-glass phase is calculated, numeri¬ 
cally exactly, as di = 2.520 for a family of hierarchical lattices, from an essentially exact (correlation 
coefficent = 0.999999) near-linear fit to 23 different diminishing fractional dimensions. To obtain 
this result, the phase transition temperature between the disordered and spin-glass phases, the cor¬ 
responding critical exponent yr, and the runaway exponent yn of the spin-glass phase are calculated 
for consecutive hierarchical lattices as dimension is lowered. 
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I. INTRODUCTION 

Singular phase diagram behavior as a function of spa¬ 
tial dimensionality d compounds the interest and chal¬ 
lenge of the phase transitions problems, as effectively 
pausing the ” phase transition of phase transitions” prob¬ 
lem. Most visible are the lower-critical dimensions, which 
are the spatial dimensional thresholds for different types 
of orderings. For example, the lower-critical threshold for 
ferromagnetic ordering in magnetic systems is = 1 for 
one-component (Ising) spins and = 2 for spins with 
more than one component. Similarly, the lower-critical 
dimensions for ferromagnetic ordering under quenched 
random fields [IH3 are respectively d/, = 2 and d^ = 4, 
for one-component spins and for spins with more than 
one component. 

A very recent experimental study Q on Ge:Mn films 
has shown the spin-glass lower-critical temperature to be 
2 < dL < 3. This is consistent with the earlier theoret¬ 
ical result of dL = 2.5 from replica symmetry-breaking 
mean-field theory. Q A numerical fit to the spin-glass crit¬ 
ical temperatures for integer dimensions has also suggest- 
ted dL = 2.5. jl^ Other theoretical work have claimed 
dL = 4 from earlier ordered-phase stability studies [n^ 
[T^ . 2 < dL < 3 from transfer-matrix studies fl3 .and 
more recently — 2 from Monte Carlo [H, Il6j| and 
ground-state studies [l^]- Renormalization-group work, 
on in effect two hierarchical lattices different from 
ours have earlier obtained 2 < dL < 3 and on a family of 
hierarchical lattices [l^, again different from ours, find 
dL close to 2.5. 

The lower-critical dimensions need not be integer, in 
view of physical systems on fractal/hierarchical lattices 
and algebraic manipulations that analytically continue. 
In fact, it would be highly interesting to find a lower- 
critical dimension that is neither an integer, nor a simple 
fraction. Our current study indicates that this is in fact 
the case for the family of hierarchical lattices studied 
here, with dL = 2.520. We obtain this result from a re¬ 
markably good fit to the renormalization-group runaway 
exponent yn from the numerically exact renormalization- 
group solution of a family of 23 hierarchical models with 
non-integer dimensions d = 2.46, 2.63, 2.77, 2.89, 3.00, 


3.10, 3.18, 3.26, 3.33, 3.40, 3.46, 3.52, 3.58, 3.63, 3.68, 
3.72, 3.77, 3.81, 3.85, 3.89, 3.93, 3.97, 4.00. Our result 
is also consistent with the results that are graphically 
displayed in Ref. [l^ for a different family of hierarchical 
lattices. 


II. LOWER-CRITICAL DIMENSION FROM 
SEQUENCED HIERARCHICAL MODELS 

Hierarchical models are constructed (20l - [^ by imbed¬ 
ding a graph into a bond, as examplified in Fig. 1, and re¬ 
peating this procedure by self-imbedding infinitely many 
times. This procedure can also be done on units with 
more than two external vertices, e.g., the layered Sier- 
pinski gasket in Ref. [1^. When interacting systems 
are placed on hierarchical lattices, their renormalization- 
group solution proceeds in the reverse direction than 
the lattice build-up just described, each eliminated el¬ 
ementary graph generating a renormalized interaction 
strength for the ensuing elementary bond. Hierarchi¬ 
cal lattices were originally introduced as present¬ 
ing exactly soluble models with renormalization-group 
recursion relations that are identical to those found in 
approximate position-space renormalization-group treat¬ 
ments of Euclidian lattices [2^ [131 , identifying the latter 
as physically realizable approximations. However, from 
the above, it is clear that any graph (or graphs 0) 
may be chosen in the self-imbedding procedure and one 
need not be faithful to any approximate renormalization- 
group solution. Hierarchical lattices have been used to 
study a varie ty of sp in-glass and other statistical mechan¬ 
ics problems. |28l460l | 

The length rescaling factor 6 in a hierarchical lattice 
is the number of bonds on the shortest distance between 
the external vertices of the elementary graph which is 
replaced by a single bond in one scale change. The vol¬ 
ume rescaling factor b'^ is the number of bonds inside the 
elementary graph. From these two rescaling factors, the 
dimensionality d is extracted, as exemplified in Fig. 1. 
In our study, 6 = 3 is used in order to treate the fer¬ 
romagnetic and antiferromagnetic correlations on equal 
footing. The lower-critical dimension of spin-glass sys- 



2 


terns is studied here by considering a systematic family 
of hierarchical lattices in all its possible decreasing di¬ 
mensions. 


III. THE SPIN-GLASS SYSTEM AND THE 
RENORMALIZATION-GROUP METHOD 

The Ising spin-glass system is defined by the Hamilto¬ 
nian 

= ( 1 ) 

<b) 


where (3 = 1/kT, at each site i of a lattice the spin 
Si = ±1, and (ij) denotes that the sum runs over all 
nearest-neighbor pairs of sites. The bond strengths Jij 
are +J > 0 (ferromagnetic) with probability 1 — p and 
—J (antiferromagnetic) with probability p. 

The renormalization-group transformation is achieved 
by a decimation, 

^ ^ ^ ^Jij SiSj-\-JSj Sf^-\-Sm 

where the additive constants Gij are unavoidably gener¬ 
ated, followed by n bond movings, 

n 

k=l 

The starting bimodal quenched probability distribu¬ 
tion of the interactions, characterized by p and described 
above, is not conserved under rescaling. The renormal¬ 
ized quenched probability distribution of the interactions 
is obtained by the convolution 


= J 


^ 3 


dJijPiJij) 


S{JL, - R{{J.,})), (4) 


where Ri{Jij}) represents the decimation and bond mov¬ 
ing given in Eqs.(2) and (3). For numerical practicality, 
the bond moving and decimation of Eqs.(2) and (3) are 
achieved by a sequence of pairwise combinations of inter¬ 
actions, each pairwise combination leading to an interme¬ 
diate probability distribution resulting from a pairwise 
convolution as in Eq.(4). The probability distribution is 
represented by 200 histograms [2^, [HI, [sj, lH, , 
which are apportioned in J ^ 0 according to total prob¬ 
ability weight. The histograms are distributed in the 
interval J+ ± 2.5(T+, where J+ and cr+ are the average 
and standard deviation of the J > 0 interactions, and 
similarly for the J < 0 interactions. 

The different thermodynamic phases of the system are 
identified by the different asymptotic renormalization- 
group flows of the quenched probability distributions. 
For all renormalization-group flows, inside the phases and 
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FIG. 1: (a) The construction of the family of hierarchical lat¬ 
tices used in this study. Each lattice is constructed by repeat¬ 
edly self-imbedding the graph. The graphs here are n parallel 
series of 6 = 3 bonds. The dimension d= 1 -f In n/ In fe of each 
lattice is given. The renormalization-group solution consists 
in implementing this process in the reverse direction, for the 
derivation of the recursion relations of the local interactions. 
The lattices shown here and 19 other lattices with the nearby 
fractional dimensions are used in our calculations, (b) The 
family of hierarchical lattices with ni parallel b — 3 series 
of 712 parallel bonds. The resulting hierarchical models are 
equivalent to the family in (a) with n = nin 2 , with respect to 
identical critical exponents and phase diagram topology in¬ 
cluding the occurrence/nonoccurrence of a spin-glass phase. 
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FIG. 2: (Color online) Critical temperatures 1/Jc and criti¬ 
cal exponents yr of the phase transitions between the spin- 
glass and paramagnetic phases as a function of dimension 
d, for the hierarchical models with antiferromagnetic bond 
concentration p = 0.5. The runaway exponents ya of the 
spin-glass phase are also shown and give a perfect fit to 
yn = -1.30908 -E 0.528513d - 0.00354805d^ leading with a 
small extrapolation to the lower-critical dimension d — 2.520 
for = 0, with a very satisfactory correlation coefficient of 
= 0.999999. 


on the phase boundaries, Eq.(4) is iterated until asymp¬ 
totic behavior is reached. Thus, we are able to calcu¬ 
late phase transition temperatures and, by linearization 
around the unstable asymptotic fixed distribution of the 
phase boundaries, critical exponents. Similar previous 
studies, on other spin-glass systems, are in Refs. [23 - 13^ . 


IV. DIMINISHING CRITICAL, RUNAWAY 
EXPONENTS, CRITICAL TEMPERATURES 
AND THE LOWER-CRITICAL DIMENSION OF 
THE SEQUENCE 

For our chosen sequence of hierarchical systems (Fig. 
1), we have calculated, at antiferromagnetic bond concen¬ 
tration p = 0.5, the phase transition temperature 1/Jc 
where the renormalization-group flows bifurcate between 
the disordered-phase and the spin-glass-phase attractor 
sinks. The spin-glass sink is characterized by an interac¬ 
tion probability distribution P(Jij) that is symmetric in 
ferromagnetism-antiferromagnetism (J^ ^ 0) and that 
diverges in interaction absolute value: The average in¬ 
teraction strength < |J| > across the system diverges 
as where n is the number of renormalization-group 
iterations and yn > 0 is the runaway exponent. The 
spin-glass sink and simultaneously the spin-glass phase 
disappears when the runaway exponent yn reaches 0.[3^ 
The calculated spin-glass phase transition temperatures, 
critical and runaway exponents are given in Fig. 2 and 
in Table I as a function of spatial dimension d. The lat¬ 
tice with d = 2.46, not having a spin-glass phase, is below 


Spatial 

Dimension 

d 

Critical 

Temperatures 

1/Jc 

Critical 

Exponents 

Vt 

Runaway 

Exponents 

VR 

2.630930 

0.519268 

0.098077 

0.058731 

2.771244 

0.747982 

0.188596 

0.129983 

2.892789 

0.890503 

0.253690 

0.191904 

3.000000 

1.001319 

0.313414 

0.246144 

3.095903 

1.091770 

0.361975 

0.294649 

3.182658 

1.168653 

0.393837 

0.338155 

3.261860 

1.237723 

0.425397 

0.377881 

3.334718 

1.298225 

0.451743 

0.414440 

3.402174 

1.354258 

0.476199 

0.448214 

3.464974 

1.404661 

0.495999 

0.479850 

3.523719 

1.452817 

0.513016 

0.509181 

3.578902 

1.496452 

0.531699 

0.536880 

3.630930 

1.538271 

0.549022 

0.563149 

3.680144 

1.577300 

0.562079 

0.587707 

3.726833 

1.613844 

0.573932 

0.610941 

3.771244 

1.649036 

0.585283 

0.633434 

3.813588 

1.682659 

0.594959 

0.654932 

3.854050 

1.714417 

0.605789 

0.675080 

3.892789 

1.745469 

0.616496 

0.693914 

3.929947 

1.774176 

0.623179 

0.712461 

3.965647 

1.802906 

0.630527 

0.730927 

4.000000 

1.829792 

0.638313 

0.747294 


TABLE I: Critical temperatures 1/Jc and critical exponents 
j/T of the phase transitions between the spin-glass and para¬ 
magnetic phases as a function of dimension d, for the hier¬ 
archical models with antiferromagnetic bond concentration 
p = 0.5. The runaway exponents yn of the spin-glass phase 
are also shown and give a perfect fit to j/_r = —1.30908 -E 
0.528513d — 0.00354805d^ , leading with a small extrapolation 
to the lower-critical dimension d = 2.520 for yn = 0, with a 
very satisfactory correlation coefficient of = 0.999999. 


the lower-critical dimension. For the 22 other consecutive 
lattices with a spin-glass phase, we have chosen to fit the 
runaway exponent values, since they gives an excellent, 
near-linear fit with 

PR = -1.30908 -E 0.528513d - 0.00354805d2, (5) 

with an amazingly satisfactory correlation coefficient of 
B? = 0.999999. This fit gives, with a small extrapolation, 
= 0 for d = 2.520. Note the near linearity, namely 
the smallness of the quadratic coefficient in Fq. (5). (In 
fact, a linear fit gives = 0 for d = 2.516, with a 
little less amazingly satisfactory correlation coefficient of 
i?2 = 0.999992.) 

Our calculated lower-critical dimension dp, where the 
spin-glass phase disappears at zero-temperature, is thus 
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seen to be cIl = 2.520, for the sequence of hierarchical 
lattices studied here. It is noteworthy that c?l is not an 
integer and not even a simple fraction, contrary previous 
examples of lower-critical dimensions (and even contrary 
to upper-critical dimensions, where mean-field behavior 
sets in) for other models. 

Another important quantity is the critical exponent 
yr = l/i^ > 0 of the phase transition between the disor¬ 
dered and spin-glass phases. This exponent is calculated 
from the scaling behavior of small deviations of the aver¬ 
age interaction strength from its fixed finite value at the 
unstable fixed distribution of the phase transition. The 
calculated critical exponents are also given in Fig. 2. 
As the spatial dimension is lowered, yx also approaches 
0. At the lower-critical dimension, yx reaches 0. The 
disordered-spin-glass phase transition disappears at dx, 
where the spin-glass phase disappears. 

V. CONCLUSION 

Our family of hierarchical lattices (Fig. 1) yields 
smooth and systematic behavior in all three quantities: 


the critical temperatures 1/Jc, the critical exponents yc^ 
and, eminently fitably, the runaway exponents y^. All 
three quantities yield the lower-critical temperature of 
dx = 2.520. It is noteworthy that dx is not an inte¬ 
ger and not even a simple fraction, contrary previous 
examples of lower-critical dimensions (and even contrary 
to upper-critical dimensions, where mean-field behavior 
sets in) for other models. 
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